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Abstract: We demonstrate the use of two high speed avalanche photo- 
^ , diodes in exploring higher order photon correlations. By employing the 

photon number resolving capability of the photodiodes the response to 
'^ , higher order photon coincidences can be measured. As an example we show 

experimentally the sensitivity to higher order correlations for three types of 

photon sources with distinct photon statistics. This higher order correlation 

rl^ ' technique could be used as a low cost and compact tool for quantifying the 

SpH, degree of correlation of photon sources employed in quantum information 

-y ' science. 
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1. Introduction 

Photon correlation functions, formalized by Glauber in the 1960s HI, have become fundamen- 
tal tools in characterizing the quantum statistical properties of light. From an experimental point 
of view, it is the class of normally ordered correlation functions ||2l corresponding to photon 
absorption as opposed to the anti-normally ordered class Q that feature as most accessible. Up 
until now, the most significant normally ordered photon correlation function is the second order 
correlation function, first measured by Hanbury-Brown and Twiss (HBT) lHISl. For example, 
the second order correlation function makes possible photon states to be categorized as clas- 
sical or quantum mechanical. Poissonian or bunched photon statistics would indicate classical 
behavior - or if sub-Poissonian statistics are observed then it would be concluded the photon 
field had a quantum character lH . 

Despite the strong interest in measuring the second order correlation function, comparatively 
little research to date has focused on measuring higher order correlations. Partly this is due 
to the complexity of the task at hand; characterizing higher order correlations of a given light 
source can require as many detectors as the order of the correlation function being measured 
itself Q . Nevertheless, from both a practical and fundamental point of view, characterization 
of higher order correlation functions is of considerable importance. In view of applications, for 
example quantum key distribution (QKD) [Sj, statistical knowledge of the source is vital for 
security. Furthermore, it has been shown for QKD using reaUstic single photon sources such as 
semiconductor quantum dots ll9][l0l[TT], the secure transmission distance can approach that of 
using a theoretical, true single photon source with knowledge of the source higher order photon 
correlation functions lfT2l . 

Traditionally, monitoring the photon field statistics by measuring higher order correlations 
is a time consuming process \T\. The time taken to acquire a large number of coincidences to 
achieve acceptable statistical bounds on the correlation function scales exponentially with the 
order of the correlation. This is usually impractical in many situations. Alternative techniques 



for measuring higher order correlations have been put forward such as two-photon detection 
lfT3l . spatial multiplexing of superconducting nanowire detectors lfT4ll and temporal multiplex- 
ing based on linear optics ITSl . While these techniques have their own individual merits, they 
suffer from other significant drawbacks including electrical cross-taUc, cryogenic operating tem- 
peratures and low coincidence rates. Moreover, to measure higher order correlations where the 
order is greater than about four, the above techniques may prove problematic due to scalability 
issues. We note that very recently a proposal to filter out cross-talk in multi-pixel detectors for 
accessing higher order correlations has been reported lfT6l . 

Measuring higher order correlations with just two detectors is a relatively unexplored area. 
One possible method is to make use of the detector photon number resolving capability. Most 
single photon detectors though can only differentiate between the zero photons and one or more 
photons 111 71 . A single or multi-photon event produces the same voltage pulse height, making 
photon number resolving impossible. On the other hand, photon number resolving (PNR) de- 
tectors produce voltage pulses in proportion to the number of photons. Hence, «-photon events 
can be precisely selected depending on the discrimination levels. In this paper, we present a 
proof of principle experiment for demonstrating sensitivity to higher order correlations based 
on high speed PNR detectors. We employ two gigahertz gated, PNR avalanche photodiodes 
(APDs) in a normal HBT setup ||4l|5l. This setup can measure the usual normally ordered and 
normalized second order correlation function g(2) defined as g(2) = ^"f""? where d^ and a 
are the usual photon creation and annihilation operators. The form of g{2) presented above is 
the most salient in terms of categorizing photon states. However, we demonstrate that due to 
the PNR capability of the detectors, two-photon, three photon and up to «-photon events can 
be selected using photon counting discriminators, where n is the PNR capability of the de- 
tector. Coupled with the joint photon detection of using both PNR detectors together, the setup 
is sensitive to higher order correlations. As we show below for a bunched photon source, the 
sensitivity to photon bunching is significantly improved over the standard g{2) measurement. 

2. Theory 

According to Glauber H], the definition of the nth order, normally ordered, normalized corre- 
lation function is given by: 

< a^"a" > 
g{n) = ,. , , (1) 

When « = 2, g{2) can be measured using the usual Hanbury-Brown Twiss experiment with 
a double single photon detector arrangement. Such a setup is depicted in Fig. [TJa) with two 
InGaAs APDs as the detectors. In this case, both detector discriminators are set so 0-photons 
are rejected while 1-photon and multi-photon events are sampled. This is readily achieved using 
most non-photon number resolving single photon detectors. As an example, if the source under 
investigation were purely thermal with a coherence time greater than the temporal resolution of 
the detector, a standard HBT experiment would yield g{2) = 2. Specifically, the measurement 
involves sampling the coincident counts at times smaller than the coherence time of the source 
and dividing by the accidental coincidences (obtained at times greater than the source coherence 
time). 

If the detectors are photon number resolving, then information on higher order correlation 
functions can be obtained by adjusting the discrimination levels. Sampling only 2-photon events 
by rejecting 0, 1 and > 3 photon events can be achieved with window discriminators with the 
window set to count only 2-photon events. We can define a quantity called the higher order 
coincidence, for example 7(4); 

r(4) = ..t7.7. ..t9... (2) 



7(4) clearly contains information regarding higher order correlation functions as we can re- 
write 7(4) in terms of the usual normalized correlation functions: 

In the example given above with the thermal source, the higher order coincidence rate would 
be higher than g{2) with 7(4) = 6. 

The above analysis can be extended to higher photon numbers. We can generalize Eq. (|2]i as: 

< ^+(ni+n2)^(ni+n2) > 

7(«i+«2) = J — ^ Tj — ^ (4) 

Here rtj indicates that the window discrimination levels of detector / is placed around the «,■ 
photon number states. Eq. (|4|i can also be re-written in terms of normalized correlation func- 
tions: 

7(«.+«2) = ^^^^ (5) 

g(ni}^g{n2) 

If the photon number resolution of both detectors is n, then the higher order coincidence can be 
measured up to 2«. Note that this arrangement circumvents the non-unity detector efficiency of 
the detectors and any inherent optical losses in the system. 

For simplicity and to illustrate the proof of principle experiment presented below, we use 
threshold discriminators rather than window discriminators. To incorporate that fact it is neces- 
sary to integrate over all photon number states greater than the lowest allowed photon number 
state. In this case Eq. Q is modified to read: 

In"'."" < fl^f"'+"2'fl("l+"2) > rj"'T72' 

^^" 1 + "2) = ^„_ "< ^.t«2fl«2 X flt«,fl«, > ^;'i^«2 (6) 

where «„„„ is the maximum photon number the detector can measure and T], is the / system 
single photon detection efficiency. The photon sources appropriate for use with this technique 
have detected mean photon fluxes much less than 2n,„ax- Note that Eq. (|6]l approximates Eq. (HJi 
when 77, is relatively small. This is due to the weighting of the higher order terms in the sums 
of Eq. (|6]l scaling as 77"'. 

3. Photon number resolution with InGaAs avalanche photodiodes 

Fig- m^) shows the layout of the experimental setup. The photon source under test is cou- 
pled into single mode fiber and photons propagate along the fiber to a 50:50 fiber beamsplitter 
Photons in each arm are detected by an InGaAs APD operating in a self-differencing mode 
at a clock rate of IGHz 1181 . The self-differencer permits cancelation of the strong APD ca- 
pacitive response leaving behind the previously obscured weak avalanche. This is achieved by 
dividing the electrical APD output into two arms, delaying one arm by an integer number of 
clock cycles before recombining the two arms using a sub tractive combiner. Fig. \V[h). After 
self-differencing, the avalanches from the detectors are amplified before being discriminated 
by ultra-fast (Ins) discriminators. The discriminated outputs are then sent to a two channel 
time-correlated photon counting card (correlator) for recording the photon arrival times. For 
measuring the detector avalanche voltage distribution, a high speed oscilloscope is employed 
directly on the amplified detector output. The APDs feature single photon detection efficiencies 
of around 15-20%. 

Fig.[2]shows the experimentally measured avalanche voltage distribution of a detector for two 
types of photon source. When excited with a 1550nm distributed feedback (DFB) type pulsed 



laser operating well above threshold (LAT, Fig.[TJc)), the avalanche voltage distribution shows 
a number of distinct peaks, black circles Fig. |2l These peaks correspond to different photon 
number, n in the incident light field. The avalanche voltage distribution can be satisfactorily 
modeled assuming the photon number distribution of the source is Poissonian with a mean 
detected photon flux, ji ^ 2.6, black line in Fig. |2l Each photon peak, « > 1 is assumed to 
be Gaussian and to reflect the statistical broadening due to avalanche noise. The widths of 
the photon peaks are scaled as -y/n relative to the 1 -photon peak width. Such a technique has 
successfully been employed to model photon number distributions for both InGaAs and Silicon 
APDs lfT9ll20l . The areas of each photon peak are proportional to the expected Poisson photon 
number distribution probabilities |19|. The main features of the LAT trace (black circles) in 
Fig.|2]are reproduced by the black line. 

The ability to detect differences in the photon statistics for dissimilar photon sources is now 
illustrated. The photon source was replaced with a multi-mode pulsed laser operating slightly 
above threshold with a Inm bandpass filter that serves to pass a single spectral mode close to a 
wavelength of 1550nm (FML in Fig.[nc)(i)). Here the average detected photon flux was similar 
to before with fi ~^ 2.8. 

This photon source gives rise to a rather different avalanche voltage distribution, (red circles. 
Fig. 121). The 0-photon peak is higher than for the laser source operating well above threshold, in- 
dicating the presence of more vacuum states. Conversely there is a suppression of the 2-photon 
and 3-photon number states. The photon source is expected to exhibit intensity fluctuations 
due to operation near to lasing threshold. Such fluctuation is manifest experimentally as pho- 
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Fig. 1. (a) Schematic of the experimental setup. The solid circles represent photons from 
the source under test. SD-APDl & SD-APD2: self-differencing avalanche photodiodes 
(APDs), discr.: electrical signal discriminator.(b) Electrical setup for measuring weak 
avalanches, (c) The three types of photon source employed and pictorial representations 
of their frequency spectra: (i) Filtered multi-mode laser operating slightly above lasing 
threshold (FML) (the mode in red box denotes allowed mode), (ii) DFB laser operating 
near lasing threshold (LNT) & (iii) DFB laser operated well above threshold (LAT). 
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Fig. 2. The experimentally measured distribution of avalanche voltages for a pulsed laser 
source operating well above threshold LAT (black circles) and for a pulsed, mode filtered 
non-DFB laser source operating close to threshold FML (red circles). Also shown are the 
theoretical distributions of avalanches for a Poissonian source (black line) and a partially 
bunched source (red line). The dashed lines correspond to the cross-over avalanche voltages 
for the experimental and theoretical avalanche distributions, as described in the text. Inset: 
Depicts the theoretical avalanche distributions for the Poissonian source (black line) and 
the partially bunched source (red line) up to an avalanche voltage of 0.75V. 



ton bunching 1211 . Bunching can be viewed as a departure from the independent statistics of a 
Poissonian source. The degree of bunching is accentuated through selection of a single mode 
by filtering ll22l . The mathematical form of the photon number distribution for this particular 
case is not known precisely, so we choose an analytic photon counting distribution function 
representing a linear superposition of a coherent state field with mean photon number /ij and a 
chaotic field with mean photon number ;U„ 1231 1241 . The following relations fix /ij & ji,, from 
experimentally measurable quantities, namely average detected photon flux, jj. = ji^ + pL„ and 
g{2) = iin {n„ + 2jLi.v)/ju2 + 1 ||23[24l. Using jx = 2.8 and g(2) = 1 .2 [|25], the resulting modeled 
avalanche distribution is shown in Fig.|2](red line). 

For both types of photon source, the theoretical avalanche distribution overestimates large 
avalanche heights. The origin of this effect is unclear at present. However, it is quite likely due 
to the quenching by the APD series resistance causing large avalanches from photon number 
states with « > 3 to be reduced in size (the linear model assumes a linear avalanche voltage 
dependence on photon number). The quenching effect can be somewhat quantified by compar- 
ing the experimental curves crossover point with the theoretical crossover point. These occur 
at 0.18 and 0.26V respectively, as indicated by the dashed lines in Fig.|2] Hence there is a 30% 
reduction in avalanche height at an avalanche voltage of 0.26V compared to the simple linear 
model prediction. 

4. Higher order correlations of the photon field 

Having shown that the detectors are sensitive to differences in the photon number distribution, 
we now present experimental results that detectors can be used to measure higher order corre- 
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Fig. 3. (a) Experimentally measured photon count histogram as a function of detector rela- 
tive delay for source FML. SD-APD 1 and SD-APD 2 discrimination levels are both set at a 
high value. Inset: Second order correlation measurement for source FML with g(2) = 1.2. 
SD-APD 1 and SD-APD 2 discrimination levels are both set at the lowest possible dis- 
crimination level, (b) Photon correlations, y of the three sources shown in Fig. [He): FML 
(squares), LNT (circles) & LAT (triangles). The arrows correspond to the mean positions 
of the photon number states assuming a avalanche voltage linear dependence on photon 
number. Inset: Detail of 7 for source LAT. Error bars are ± 1 standard deviation of the 
plotted (mean) values and are plotted for all data; however only some error bars are visible 
due to the size of the data points. 



lations of the photon field. The setup illustrated in Fig.lTJa) is used. We stress that the setup will 
be exactly the same as a conventional one that measures the second order correlation function 
g(2) if the photon number resolving function is disabled. By setting the discrimination level at 
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Fig. 4. Simulated correlation value as a function of SD-APD 2 discrimination level for SD- 
APD 1 discriminator level set at ;? = 7 photon number state. FML (squares), LNT (circles) 
& LAT (triangles). 



300m V, a value far greater than single photon generated avalanches, for both detectors SD-APD 
1 and SD-APD 2, the measurement gives a correlation result that is different from g(2). Plotted 
in Fig.[3la), is a histogram for source FML when both the discrimination levels of SD-APD 1 
and SD-APD 2 are at high values. At zero time delay (Af = 0) a prominent correlation peak is 
observed due to photon bunching. The correlation here is given by Eq. (4), y. The correlation 
value is measured as the ratio of peak height at t=0 to the average height at Af 7^ 0. 7 is much 
higher than the the equivalent g{2) value of g(2) — 1 .2 , as shown by the histogram in the inset 
of Fig. la). 

For different discrimination levels, correlation plots similar to Fig. 13 a) are collected and the 
correlations, 7 evaluated. The detector discriminator level of detector SD-APD 1 was kept at 
300m V. Fig.[3lb) plots the resulting correlations, 7 as a function of SD-APD 2 discrimination 
level for the three sources. All three sources show 7 > 1 even at the lowest discrimination level 
used. This can be readily understood in terms of elevated photon bunching through higher or- 
der correlations. The «''' order correlation is the expectation of the joint detections of n photons 
correlations. For photons that are indistinguishable and statistically dependent, there is a facto- 
rial increase of the available permutations of photon amplitudes as n rises Il26ll27l . At higher 
discrimination levels, higher photon number states feature more markedly than lower photon 
number states and the correlation is therefore expected to be higher than for lower photon num- 
ber states. This is in direct contrast to what would be expected for the case of statistically inde- 
pendent photons. In this case the photon statistics are Poissonian and no bunching is possible 
for any order of correlation. 

Photon source LAT shows the least amount of correlation; this is expected for a laser source 
operating well above threshold with largely Poissonian statistics [2|. The correlation for LAT is 
relatively flat over the entire discrimination range of SD-APD 2. 7 for this source attains ~ 1 .3 
at the very highest D2 discrimination level employed, Fig.[3jb), inset. Operating the same laser 
near to threshold (LNT) shows a marked increase in correlation which is due to increased 
intensity fluctuations leading to increased photon bunching. Source FML displays the highest 
correlation values; attaining 7^ 47.1 ±4.5 for these measurements. Intensity fluctuations for 



this source are by far the highest of all sources under test. 

All three sources under study feature different photon number statistics which is incorpo- 
rated in our model through the mean photon numbers /ij & ji„. These mean photon numbers 
are derived from the total average photon number ji and g{2). The second order correlation 
function g{2) was measured separately for each of the sources FML, LNT & LAT with the 
values g{2) = 1.2, 1.075 & 1.001 respectively (the inset of Fig. 13 a) shows the measured g(2) 
photon coincidence histogram for source FML). 

As remarked previously, avalanche height self-limiting due to the APD series resistance is 
likely to distort the linear avalanche voltage dependence on photon number for high photon 
number. We can incorporate this fact by increasing n,„„t in the model for the detectors. Fig. 
|4] shows the simulated results based on Eq. (|6]l for y plotted with up to «„,„, = 7 photon de- 
tection events . Each point corresponds to selecting exactly n avalanche detections based on a 
linear dependence of the avalanche peaks (as shown by the arrows in Fig.[3jb)). The x-axis is 
plotted reciprocally to emphasize the avalanche self-limiting. 7 for all sources shows a similar 
growth and dependence to that observed experimentally, qualitatively confirming our underly- 
ing model. 

Finally, we remark that although the agreement between the experimental and theoretical 
data is good, it is only qualitative. The photon number resolving capability of the detector 
worsens at higher photon number due to avalanche broadening. Avalanche broadening causes 
adjacent photon number avalanches to overlap. This fact is not incorporated in the model for 
simplicity. The next step is to realize a SD-APD with narrow photon number avalanche peaks 
having minimal overlap; similar to the performance of visible light photon counters (VLPCs), 
for example 1281 . This will permit efficient window discrimination between n and « ± 1 photon 
number avalanche peaks. Then the correspondence between the experimentally measured 7 and 
Eq. (4) would be direct. 

5. Conclusion 

In summary we have demonstrated a unique, proof of principle technique for inferring higher 
order photon correlations based on high speed (GHz) photon number resolving APDs. Using 
three light sources with differing photon statistics we have shown that higher order correlations 
of these sources can be detected. Moreover, accessing these higher order correlations illustrates 
the sensitivity of the method when trying to distinguish between sources compared to a usual 
g{2) measurement. We believe the technique will be of use and importance in characterizing 
photon sources involved in the rapidly evolving field of quantum information science. 
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